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We put forward a scheme for controlling Josephson-like tunneling in two- component Bose-Einstein 
condensates coupled with microwave field via Feshbach resonance and tuning aspect ratio of trapping 
potential. We prove how to realize a perfect periodic oscillation from a fast damped and irregular 
oscillation on relative number of atoms in future experiment. In particular, intensity of Josephson- 
like tunneling can be successfully controlled through controlling speed of recovering the initial value 
of intra-atomic interaction and aspect ratio of trapping potential. Interestingly, we find that relative 
number of atoms represents two different types of oscillation in respond to periodic modulation of 
attractive intra-atomic interaction. 
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The existence of a Josephson current is a direct man- 
ifestation of macroscopic quantum phase coherence [l| 
and have numerous important applications in condensed 
matter physics, quantum optics and cold atom physics, 
for example, precision measurement, quantum compu- 
tation. The physical origin of the Josephson current is 
the temporal interference of the two systems which must 
both have a well defined quantum phase and a differ- 
ent average energy per particle, respectively. Recently, 
the experimental realization of multi-component Bose- 
Einstein condensates (BECs) of weakly interacting alkali 
atoms hasprovided a route to study Josephson effect 
@, 0, H S IE 0, H, i, in a controlled and tunable way 
by means of Feshbach resonance [ll|, [l2|, [l3|, \i% [l5|, LL6[ so 
far unattainable in charged systems. The physical origin 
of the low-energy Feshbach resonances is the low-energy 
binary collisions described by the difference of the initial 
and intermediate state energies which can be effectively 
altered through variations of the strength of an external 
magnetic field. 

In this Letter, we put forward a scheme for controlling 
Josephson-like tunneling in two-component BECs cou- 
pled with microwave field via Feshbach resonance and 
tuning aspect ratio of trapping potential, we find that, 
through time-dependent tuning attractive intra-atomic 
interaction, a perfectly periodic oscillation can be ob- 
tained successfully. Especially, by controlling speed of 
recovering the initial value of attractive intra-atomic in- 
teraction and aspect ratio of trapping potential, we can 
successfully control intensity of Josephson-like tunneling. 
Furthermore, relative number of atoms represents two 
different types of oscillation in respond to periodic mod- 
ulation of attractive intra-atomic interaction. 

We consider a system of longitudinal elongated two- 
component BECs coupled by the microwave field with 
the effective Rabi frequency Vt and finite detune 5, where 
the Q and 5 are independent upon time and space coordi- 



nate as experimental case [3[. At zero temperature, this 
system can be described by macroscopic wave function ipj 
(normalized to unity, i.e., J(\^i\ 2 + |^2| 2 )^ = 1) which 
obey the dimensionless one-dimensional Gross-Pitaevskii 
equations (GPEs) as [a,uj 

2 

[G jj (t)\ip j \ 2 + G jk (t)\ip k \ 2 ]ip j (z,t) + 

( _i);+iA^ (M) + _% fc(M)) (i) 

where j, k = 1,2, j ^ k, and co± is radial harmonic trap- 
ping frequency, time and coordinate unit are 2/co_\_ and 
a±_ = y/h/muj_\_, respectively, and V(z) is the longitudi- 
nal confining potential, Gjj = (ANhdjj) / \ma\uj±) and 
Gjk = (4:Nhdj k)/ (jna are dimensionless intra- and 
inter- species effective atomic interaction potentials, re- 
spectively, where ajj and ajk are corresponding s-wave 
scattering lengths, respectively, N is the total atomic 
number, and m is the atomic mass assumed the same 
for both two species. 

In the following discussion, we consider two-hyperfine 
spin states of 7 Li bose atom which the magnitude and 
sign of the atomic interaction can be tuned to any value 
via a magnetic- field dependent Feshbach resonance. And 
we choose dimensionless V(z) = 7 2 z 2 , where 7 = uj z /uj±_ 
is aspect ratio of trapping potential, uo z is longitudinal 
harmonic trapping frequency and oj±_ = 2tt x 625 Hz (so, 
the time and coordinate units are 0.51 ms and a± = 1.51 
fim respectively). Moreover, considering the condensates 
will be unstable and collapse when the number of parti- 
cles is large enough [17j, so we choose N = 6000, which 
provides a safe range of parameters for avoiding instabil- 
ity and collapse occurring. And the effective Rabi fre- 
quency and finite detune are chosen as Q = l.buu± and 
S = 0.5cu_\_ considering experimental feasibility. 
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FIG. 1: (Color online) Temporal evolution of atomic density 
and relative number of atoms Nd for initial Gaussian wave 
function in response to various aspect ratio of trapping po- 
tential 7. (a) Corresponding to attractive atomic interactions 
for the parameters as Gu = G22 = 4Ci2 = —10. (b) Corre- 
sponding to repulsive atomic interactions for the parameters 
as Gu = G22 = 4<3i2 = 10. 



First of all, we consider the simplest case, i.e., intra- 
and inter-species effective atomic interaction potentials 
are time-independent attractive or repulsive case, which 
can be easily realized using Feshbach resonance experi- 
mentally. Here, we have chosen an = <222 is just con- 
sidering the experimental feasibility and assumed that 
tuning atomic interaction wouldn't destroy the conden- 
sates in this system. We denote relative number of atoms 
between two species as A^, i.e., Nd = Ni — 7V 2 , where 
N± = f\il)i\ 2 dz and N 2 = J \^2\ 2 dz, respectively. We 
study Josephson-like tunneling in this case and the re- 
sults are shown in Fig. 1, which is based on numerical 
simulation of Eqs. (pQ) for the initial Gaussian wave func- 
tion under the various aspect ratio of trapping potential. 
As shown in Fig. la corresponding to attractive atomic 
interactions, during the evolution process, exchange of 
atoms is always damped very fast and almost suppressed 
completely irrespective of aspect ratio of trapping poten- 
tial. But, on the contrary, in repulsive atomic interac- 
tions corresponding to Fig. lb, exchange of atoms vary 
significantly in response to various aspect ratio of trap- 
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FIG. 2: (Color online) Temporal evolution of Nd for initial 
Gaussian wave function in response to various aspect ratio of 
trapping potential 7, which is based on numerical simulation 
of Eqs. © for the case E (77) =0. (a) Corresponding to 
attractive atomic interactions for the parameters as Gu — 
G22 = 4Gi2 = —10. (b) Corresponding to repulsive atomic 
interactions for the parameters as Gu = G22 = 4Ci2 = 10. 

ping potential. Detailed analyses of above phenomena 
will be given in the following. 

From Eqs. (pQ) and taking boundary condition into 
account, it is not difficult to derive following equations, 

dN d 4n r 

—j- = — / WP1P2 sin 4>d)dz, (2) 
at uj±_ J 

where wave functions have been chosen as ipi = 
v /pTexp(z0i), ^ 2 = v /p 2 "exp(^ 2 ), respectively, and <j> d = 
4>i~ 4>2- From Fig. la, we can see that maximum pi will 
be much smaller than maximum p2 very quickly for all 
aspect ratio of trapping potential, so maximum ^pip2 
will decrease rapidly. Moreover, if considering well local- 
ized density profiles of two species as shown in Fig. la, 
from Eqs. ([2j), we can immediately conclude that am- 
plitude of oscillation of Nd will be damped very fast and 
suppressed completely in this case. But from Fig. lb cor- 
responding to repulsive atomic interactions, we can see 
that the atomic densities of two species vary significantly 
in response to various aspect ratio of trapping potential 
and represent nonlocalized profiles, so time evolution of 
amplitudes of N d represents very different behavior from 
attractive case. 

As pointed out above, Nd will represent a very irregular 
and fast damped oscillation in time-independent attrac- 
tive atomic interactions and represent a more irregular 
oscillation in time-independent repulsive case. But, from 
the experimental point of view, valuable and applicable 
systems usually have regular and stable characteristics 
rather than irregular and unstable one. Fortunately, we 
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FIG. 3: (Color online) Temporal evolution of relative number 
of atoms Nd for initial Gaussian wave function in response to 
various aspect ratio of trapping potential 7. (a) Correspond- 
ing to attractive atomic interactions for the parameters as 
G12 = -2.5, r = 105, ti = 0.5t 2 = 280. (b) Corresponding to 
repulsive atomic interactions for the parameters as G12 =2.5, 
r = 105, ti = 0.5t 2 = 280. 

find out that, through time-dependent tuning atomic in- 
teraction via Feshbach resonance, regular and stable os- 
cillation of Nd can be successfully realized in the attrac- 
tive case, but not be realized in the repulsive case. Now, 
let us to understand physics from both theoretical anal- 
ysis and numerical simulation aspects. 

In order to study how to control Josephson-like tun- 
neling in this system, we diagonalize linear parts of Eqs. 
(PQ) by means of applying an unitary transformation U 
on both sides of Eqs. (pp), then it is transformed into 
following form, 

ig- t <P±(z,t) = [—^ + V(z)]<p±(z,t)+(t\<p±\ 2 <p±(z,t) 
+C2 ± I^tIV±(M) + (3) 

where 

tj ( V>i W v+e^ \ = (*+\ TT= ( s ™° ~ cos0 \ 
\1p2J V <P-e~ ivt )~ \ *- J \cos0 sin6» y 

here (f = [(3Gn + 2G i2 + 3G 22 ) T 4(Gn - G 22 ) cos 29 + 
(Gn - 2Gi2 + G 22 )cos40]/8, g = [(G n + 2G i2 + 
G 22 ) - (G n - 2G 12 + G 22 )cos40]/4, C 3 ± = [2(G n - 
G 22 )sin2^T(Gii-2G 12 + G 22 )sin4^]/8 ) C 4 ± = [2(G n - 





'■ : < ' ", : ;'■ ,'' 'i/,, 11 


I 

t 2 = 112 _ 
t 2 =84 . 

I 11 

- - $«Z 






40 


f 80 (3) 120 160 


ii , i i , i 


=1/250 



40 t 80 (b) 120 160 



FIG. 4: (Color online) Temporal evolution of Nd for initial 
Gaussian wave function in response to (a) various speed of 
recovering the initial value of intra-atomic interaction for the 
parameters as 7 = 0, G12 — —2.5, r = 21, t\ — 56, and (b) 
various aspect ratio of trapping potential for the parameters 
as G12 = -2.5, r = 21, ti = 0.83*2 = 56. 

G 22 )sin2#T(Gn -2Gi 2 + G 22 )sin4#]/4, C 5 ± = [(Gn - 
2Gi 2 + G 22 )sin 2 2#]/4, C 6 ± = [2(Gn - G 22 )sin2# ± 
(Gn - 2Gi 2 + G 22 ) sin40]/8, and = (1/2) tan" 1 (Q/S), 
77 = [(5/uj_) 2 + (Q/uj±) 2 ] 1 / 2 . Now, we can see that, com- 
paring with the original Eqs. (pp), we introduce a new 
term denoted by £(77) which is exponential time depen- 
dence, and it is not included in regular coupling nonlin- 
ear Schrodinger equations. It must be pointed out that 
11(77) depends not only on the effective atomic interac- 
tion potentials Gjk (j,k = 1,2), but on the effective Rabi 
frequency ft and finite detune <S, i.e., 6 and 77. 

From the expression of £(77), we can find that, if 
tp±(z,i) vary slowly over time period 71-/77, then the ef- 
fects of 11(77) on Josephson-like tunneling can be averaged 
out over time period 71-/77 during the evolution process. 
But if ip±(z,i) vary significantly over time period 71-/77, 
then £(77) will strongly affect Josephson-like tunneling. 
In the following, we will investigate effects of £(77) on 
Josephson-like tunneling in detail. 

Firstly, we consider the special case that £(77) = 0. 
From Eqs. ([3]) and transformation relation between orig- 
inal and transformed wave function, we can obtain that 
Nd will represent a large-amplitude periodic oscillation in 
both attractive and repulsive atomic interactions. Above 
conclusions are well confirmed by Fig. 2, which is based 
on numerical simulation of Eqs. (|3]) and then applying 
inverse unitary transformation £7 _1 on \£± in order to 
obtain original wave function. So, it is reasonable to 
infer that the reasons for fast damped and irregular os- 
cillation of Nd are mainly due to the effects of £(77) on 
Josephson-like tunneling. 

But, key question is how to efficiently suppress effects 
of £(77) on Josephson-like tunneling in order to obtain 
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perfectly periodic oscillation of Nd- In the following, we 
prove that, only in the attractive case, through time- 
dependent tuning intra-atomic interaction via Feshbach 
resonance, we can successfully control the effects of £(77) 
and obtain a perfectly periodic oscillation on relative 
number of atoms. Considering the experimental feasi- 
bility, in the following case, we choose time-dependent 



intra-atomic interaction as 



Gi 2 exp(t/r), 
Gjj(t) = { G 12 exp[(t 2 -£)/r], 



h<t ^ t 2 - 
t > t 2 . 



(4) 



where j =(1, 2), Gn(t) = G 22 (t) and G\ 2 is time- 
independent parameter determining inter-atomic inter- 
action, r determines speed of tuning intra-atomic inter- 
action depending on the values of t\ and t 2 at fixed min- 
imum Gjj(t). The results are shown in Fig. 3. It is very 
interesting that Nd finally represents a perfectly periodic 
oscillation in the attractive case, but not in the repulsive 
case. Above phenomena are mainly due to successfully 
suppress the effects of £(77) on Josephson-like tunneling 
in the attractive case, but not in the repulsive case, which 
underlying physics is just because (p±(z,t) vary signifi- 
cantly over time period tt/t] in the repulsive case. So, 
in the following, we will concentrate on attractive case. 
It is important to point out that, in the attractive case, 
the amplitude of Nd is strongly dependent upon aspect 
ratio of trapping potential as shown in Fig. 3a. Based 
upon the above analysis, we put forward a scheme for 
controlling the amplitude of oscillation through control- 
ling speed of recovering the initial value of intra-atomic 
interaction (i.e., changing value of t 2 but fixing other pa- 
rameters in Eqs. (|4|)) and tuning aspect ratio of trapping 
potential. 

The results are shown in Fig. 4. It is very interesting 
to note that oscillatory amplitude of Nd can be controlled 
successfully not only via Feshbach resonance but tuning 
aspect ratio of trapping potential. Above phenomena can 
be well understood as follows: If ip±(z, t) vary slowly over 
time period 71-/77, then the effects of £(77) on Josephson- 
like tunneling can be effectively suppressed, and Nd will 
represent a perfectly periodic oscillation which amplitude 
is determined by atomic density of two species. 

We also study Josephson-like tunneling in response 
to periodic modulation of attractive interaction [l2[ 
and various aspect ratio of trapping potential. In this 
case, attractive atomic interactions have been chosen 
as Gu(t) = G 22 (t) = G 12 [l + 8sin 2 (0.01t)], where 
G12 = —2.5. As shown in Fig. 5, it is very interest- 
ing to note that Nd represents two different types of os- 
cillation, which one represents small frequency behavior 
originated from periodic modulation of attractive inter- 
action, but the other represents large frequency behavior 
originated from coupling microwave. Above effects pro- 
vide a promising protocol to realize a new type of switch 
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FIG. 5: (Color online) Temporal evolution of relative number 
of atoms Nd for initial Gaussian wave function in response to 
periodic modulation of attractive intra-atomic interaction and 
various aspect ratio of trapping potential for the parameter 
as G12 — —2.5. 



controlled by periodically managed Feshbach resonance 
in future experiments. 

In conclusion, we prove how to successfully control 
Josephson-like tunneling in attractive two-component 
BECs coupled by microwave field via Feshbach resonance 
in respond to various aspect ratio of trapping potential. 
Recent developments of controlling the scattering length 
and realization of multi-component BECs in the exper- 
iments allow for the experimental investigation of our 
prediction in future. 
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